Modified Spin- Wave Theory of Nuclear Magnetic Relaxation in One-Dimensional 
Quantum Ferrimagnets: Three-Magnon versus Raman Processes 
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Nuclear spin-lattice relaxation in one-dimensional Heisenberg ferrimagnets is studied by means of 
a modified spin-wave theory. Calculating beyond the first-order mechanism, where a nuclear spin 
directly interacts with spin waves through the hyperfine coupling, we demonstrate that the exchange- 
scattering-enhanced three-magnon nuclear relaxation may generally predominate over the Raman 
one with increasing temperature and decreasing field. Recent proton spin-lattice relaxation-time 
(Ti) measurements on the ferrimagnetic chain compound NiCu(C7H6N20s)(H20)3-2H20 suggest 
that the major contribution to 1/Ti be made by the three-magnon scattering. 

PACS numbers: 75.10.Jm, 75.50.Gg, 76.50.+g 



INTRODUCTION 



Nuclear magnetic resonance (NMR) is an effective 
probe to the collective motions of electronic spins and 
therefore we take a great interest in microscopically in- 
terpreting it. The spin- wave formalism has played a cru- 
cial role in this context. Van Kranendonk, Bloomi and 
Moriya^ made their pioneering attempts to describe the 
nuclear spin-lattice relaxation time T\ in terms of spin 
waves. Oguchi and Keffer— further developed the spin- 
wave analysis considering the three-magnon nuclear re- 
laxation mechanism as well as the Raman one. Pincus 
and Beeman^& claimed that the exchange correlation be- 
tween spin waves should significantly accelerate the nu- 
clear spin relaxation. The spin-wave excitation energy 
is usually much larger than the nuclear resonance fre- 
quency and thus the single-magnon relaxation process is 
rarely of significance. The Raman process consequently 
plays a leading role in the nuclear spin-lattice relaxation. 
Because of the (4S) -damping factor to the Holstein- 
Primakoff magnon series expansion, the multi-magnon 
scattering is much less contributive within the first-order 
process, where a nuclear spin directly interacts with spin 
waves through the hyperfine coupling. However, the 
second-order process, where a nuclear spin flip induces 
virtual spin waves which are then scattered thermally 
via the four-magnon exchange interaction, may generally 
enhance the relaxation rate. This is the fascinating sce- 
nario written by Pincus and Beeman. 

It is unfortunate that their spin-wave nuclear relax- 
ation theory is not effective in low dimensions but is 
valid far below the transition temperature. The con- 
ventional spin-wave theory applied to low-dimensional 
magnets ends in failure with diverging magnetizations. 
In such circumstances, Takahashii gave a fine descrip- 
tion of the low-dimensional ferromagnetic thermody- 
namics in terms of modified spin waves. His idea of 
introducing a constraint on the magnetization so as 
to control the number of spin waves was further ap- 
plied to antiferromagnets&2iiii and ferrimagnets>ii*i2ii^ 
Even frustrated antiferromagnets^ii£*A£ and random- 



bond ferromagnetsii were discussed within this renewed 
spin- wave scheme. 

The ferrimagnetic modified spin-wave theory is par- 
ticularly useful in illuminating both statioi2*i2*i^ and 
dynamio22i2ii2£ properties. One-dimensional ferrimag- 
nets have lately attracted much attention especially in 
the context of designing molecule-based ferromagnets. 
Assembling molecular bricks in such a way as to obtain 
a low-dimensional system with a nonzero resultant spin 
in the ground state and then coupling the chains or the 
layers again in a ferromagnetic fashion, one can in prin- 
ciple obtain a molecular magnet. A series of bimetallic 
chain compounds 2 ^' were synthesized in such a strategy. 
Another approaches' to molecular magnets consists of 
bringing into interaction metal ions and stable organic 
radicals. Homometallic materials such as tetrameric 
bond-alternating chain compounds^ and trimeric inter- 
twining double-chain compounds— are distinct ferrimag- 
nets of topological origin. Such synthetic endeavors have 
stimulated several experimentalists 29 ' 30-31 to measure T\ 
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FIG. 1: Modified spin-wave calculations of the specific heat 
C and the magnetic susceptibility \ as functions of temper- 
ature for the ferrimagnetic Heisenberg chains. The original 
(Takahashi) and our (Yamamoto) schemes are compared with 
quantum Monte Carlo (QMC) calculations. 
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FIG. 2: Diagrammatic representation of various nuclear spin-lattice relaxation processes. Solid arrows, designating spin waves 
which are emitted in the first-order processes, induce a nuclear spin flip (x) via the hyperfine interaction, while broken arrows, 
depicting the four-magnon exchange correlations, thermally scatter the first-order spin waves as virtual excitations, where spin 
waves of ferromagnetic and antiferromagnetic aspect are distinguishably drawn by straight and wavy arrows, respectively, (a) 
The first-order direct (single-magnon) relaxation processes; (b) The first-order Raman (two-magnon) relaxation processes; (c) 
The first-order and second-order three-magnon relaxation processes, where q = — k^, are related to each other through nonlinear 
equations and are therefore inseparable. Considering the nuclear-electronic energy conservation, processes in solid and dotted 
frames are of great and little significance, respectively, whereas those in broken frames are relevant according to the constituent 
spins S and s. Labels W <TtT and W" 7 " on feasible processes are explained in Appendix ITU 



on ferrimagnetic chain compounds. Thus motivated, 
we have started a modified spin-wave exploration^ of 
the nuclear spin dynamics in one-dimensional Heisen- 
berg ferrimagnets beyond the conventional Raman mech- 
anism. Here we give a full description of the the- 
ory and finally show a strong evidence of the proton 
spin relaxation in the ferrimagnetic chain compound 
NiCu(C 7 H e N20 e )(H 2 0) 3 ■ 2H 2 being mediated by the 
three-magnon scattering rather than the Raman one. 

II. MODIFIED SPIN- WAVE SCHEME 

We consider ferrimagnetic Heisenberg chains of alter- 
nating spins S and s, as described by the Hamiltonian 

N 

n = J2[ JS n-(sn-i + s n )-(g s S z n + g s s z n ) f i B H]. (2.1) 

71=1 

Introducing bosonic operators for the spin deviation in 
each sublattice via = {2S — a\a i ) 1 / 2 a i , S z — S — a\ai, 
sf = b\(2s — fej^i) 1 / 2 , s* = — s + b\bi, and assuming that 
O(S) = O(s), we expand the Hamiltonian with respect 



to 1/5 as U = -2SsJN + Hi + Hq + 0(5" 1 ), where 
Tii contains the 0(S l ) terms which are explicitly given 
in Appendix 1X1 



Our scheme^* of modifying the spin- wave theory is dis- 
tinct from the original idea proposed by TakahashS and 
Hirsch et al.£ In their way of suppressing the divergence 
of the sublattice magnetizations, an effective Hamilto- 
nian with a Lagrange multiplier included, instead of the 
original Hamiltonian, is diagonalized subject to zero stag- 
gered magnetization. On the other hand, we first diag- 
onalize the Hamiltonian keeping the dispersion relations 
free from temperature and then introduce a Lagrange 
multiplier in order to minimize the free energy subject to 
zero staggered magnetization. The two approaches are 
compared in Fig. ^ Our new scheme is much better at 
describing the antiferromagnetically peaked specific heat 
and the low-temperature diverging susceptibility. Our 
approach converges into the paramagnetic behavior at 
high temperatures for both the specific heat and suscep- 
tibility. 
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III. NUCLEAR SPIN-LATTICE RELAXATION 

The hyperfine interaction is generally expressed as 

H h f = 9S^B^lNl + Y,n(h A n S n + A n S n) 

+9simh 1 Kl + Tn{\B-s- + B z n s z n ), (3.1) 



where A° (B^) is the dipolar coupling tensor between the 
nuclear and nth larger (smaller) electronic spins. Since 
H.0 and Hhf are both much smaller than H. i , they act as 
perturbative interactions to the linear spin-wave system. 
If we consider up to the second-order perturbation with 
respect to V = Ho + Hhf, the probability of a nuclear 
spin being scattered from the state of I z = m to that of 
I z — to + 1 is given by 



/ V 



E 



V|to)(to|V 



Ei — E,„ 



m - E t ), 

(3.2) 

where i and / designate the initial and final states of 
the unperturbed electronic-nuclear spin system, whose 
energies are Ei and Ej, respectively. The nuclear spin- 
lattice relaxation time is then given by Ti = (I — m) (J + 
m + l)/2W. 

Equation l|3.2[l contains various scattering processes, 
which are diagrammatically shown in Fig. [21 Due to 
the considerable difference between the nuclear and elec- 
tronic energy scales, huj-^ <C J, the direct process, involv- 
ing a single spin wave, is rarely of significance. Consid- 
ering further that the antiferromagnetic spin waves are 
higher in energy than the ferromagnetic ones, uj^ < uj^ , 
at moderate fields, the intraband spin-wave scattering 

(2) 

dominates the Raman relaxation rate 1/T| , whereas 
both the intraband and intcrband spin-wave scatterings 
contribute to the three-magnon relaxation rate 1/T| . 
Within the first-order mechanism, 1/Ti is much smaller 

than 1/T; 4 However, the first-order relaxation rate 
may be enhanced through the second-order mechanism. 
We consider the leading second-order relaxation, that 
is, the exchange-scattering-induced three-magnon pro- 
cesses, as well as the first-order ones. The second-order 
Raman processes, containing two virtual magnons, are 
much more accidental due to the momentum conser- 
vation and much less contributive due to the (4S') _1 - 
damping factor in the Holstein-Primakoff magnon se- 
ries expansion. As for the four-magnon scattering, the 
first-order processes are nonexistent, whereas the second- 
order ones originate in the six-magnon exchange interac- 
tion and therefore contain two virtual magnons. Thus 
and thus, all other higher-order processes do not signifi- 
cantly change the relaxation scenario. We explicitly for- 
mulate 1/T X (2) and l/T^ in Appendix |U 



IV. THREE-MAGNON VERSUS RAMAN 
PROCESSES 

We calculate the cases of (S, s) = (1,5) an d 
(S, s) = (§,5), which are relevant to several major 
materialsi22i2£i2& Figure shows 1/Ti as a function 
of temperature and an applied field. The exchange- 
scattering-enhanced three-magnon relaxation rate gener- 
ally grows into a major contribution to 1 /Ti with increas- 
ing temperature and decreasing field. As temperature in- 
creases, n^T decreases at k ~ but otherwise increases^ 
In one dimension, excitations at k ~ predominate in 

the Raman processes, while all the excitations are effec- 

(2) (3) 
tive in the three-magnon processes. l/Tf and 1/T{ 

are hence decreasing and increasing functions of temper- 
ature, respectively, unless temperature is so high as to 
activate the antiferromagnetic spin waves. The field de- 



pendences of 1/T 1 (2; and 1/T^ are also in striking con- 
trast. At moderately low temperatures and weak fields, 
(2) 

hujjq <C k B T <C J, T{ is approximately evaluated as 
1 2[(g s A z S -g s B z s)^ B h 7N ] 2 



(3) 



(2) 



xexp 



nhSs(S — s)J 



2k B T 



(4.D 



\2k B T) 



where K is the modified Bessel function of the sec- 
ond kind and behaves as K (huj N /2k B T) ~ 0.80908 - 

(2) 

\n(huJN / k B T) . Thus the field dependence of 1/T| ; is 
initially logarithmic and then turns exponential with in- 
creasing field. Equation 1|B2J| is much less analytical 
but suggests much stronger power-law diverging behavior 
with decreasing field. Therefore, the three-magnon relax- 
ation predominates over the Raman one at weak fields. 

In Fig. 0] we plot the crossover points on which 
1/T{ ' — 1/Ti • A Raman-to-three-magnon crossover 
may generally be detected with increasing tempera- 
ture and decreasing field. The ferrimagnetic nuclear 
spin-lattice relaxation is sensitive to another adjustable 
parameter A T /B T , that is, the location of the probe 
nuclei^ At the special location of A T / B T ~ (d s /d s ) 3 ~ 
(S/s) a , where ds (d s ) is the distance between the nu- 
clear and larger (smaller) electronic spins, the a excita- 
tion mode is almost invisible to the nuclear spin^ In 
the case of (S, s) = (1, 4), for exmaple, the nuclear spin 
located as A T /B T ~ 1/2 hardly relaxes. Any T\ mea- 
surements should be performed away from such magic 
points. 

V. INTERPRETATION OF EXPERIMENTS 

We are further excited to compare our theory 
with recent experimental findings. Fujiwara and 
Hagiwara 31 measured T\ for proton nuclei in the bimetal- 
lic chain compound NiCu(pba)(H20)3.2H20 (pba = 
1, 3-propylenebis(oxamato) = C7H6N2O6) 36 comprising 
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FIG. 3: Modified spin-wave calculations of typical temperature (the left four) and field (the right four) dependences of the 
nuclear spin-lattice relaxation rate, where gs = g s = g, A T /B T = and (B~ /B z ) 2 = 4. 1/Tj 2 ' and are plotted by 

dotted and broken lines, respectively, while + 1/T{ 3 ' = l/Ti, which is observable, by solid lines. 



ferrimagnetic chains with alternating octahedral Ni 2+ 
and square-pyramidal Cu 2+ ions bridged by oxamato 
groups. The measured susceptibility 37 is well repro- 
duced with g s = 2.22, g s = 2.09 and J/k B ~ 121 K. 
Since the protons relevant to the T\ findings are located 
in close vicinity to Cu ions, we may set the coupling 
constants for A T /B T = 0. We further assume that 
B z = 1.37 x 10 20 T 2 /J and (B-/B z ) 2 = 5, which can 
be consistent with the crystalline structure. 36 

The thus-calculated l/Ti is compared with the ob- 
servations in Fig. [SJ Considering that there may be 
larger uncertainty in the experimental analysis at lower 
temperatures^ the theoretical and experimental find- 
ings are in good agreement and the slight discrepancy 
between them may be attributable, for instance, to weak 
momentum dependence of and the protons of wide 
distribution. Figure [S] (b) shows that the increasing be- 
havior of 1 /T\ at high temperatures originates from the 
three-magnon contribution. Figure Efc) more impres- 
sively demonstrates the relevance of the three-magnon 
scattering to the proton spin relaxation. The strong field 
dependence can never be explained by the Raman process. 
Since 1 /T| within the first-order mechanism stays much 
smaller than the observations, the exchange- scattering- 
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FIG. 4: The crossover point as a function of temperature and 
an applied field, where gs = g s and [B~ / B z ) 2 — 4. 



induced three-magnon process is essential in interpret- 
ing such accelerated relaxation. We are eager to have 
reliable observations at lower temperatures and weaker 
fields. More extensive NMR measurements on the re- 
lated compounds are encouraged. 



VI. SUMMARY 



There exist pioneering T\ measurements on the lay- 
ered ferromagnet CrClg^ and the coupled-chain anti- 
ferromagnet CsMnCla^HsOf 3 ^ which give evidence of 
the relevant three-magnon scattering. However, they are 
both, in some sense, classical findings under the existing 
three-dimensional long-range order. No author has ex- 
plored one-dimensional quantum ferrimagnetic dynamics 
with particular interest in multi-magnon scattering be- 
yond the Raman mechanism. We have reported the first 
evidence of the three-magnon scattering dominating one- 
dimensional nuclear spin relaxation. We hope the present 
research will stimulate further measurements and lead to 
close collaboration between theoretical and experimental 
investigations. 
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FIG. 5: Proton spin relaxation-time measurements on NiCu(pba)(H20)3-2H20 (Ref. I31T) compared with our theory, (a) 1/Ti 
as a function of temperature at various values of an applied field; (b) 1/Ti as a function of temperature at 3.15 T; (c) 1/Ti as 



a function of an applied field at 280 K. In (b) and (c), the Raman (l/T^ 2 ') 
plotted by dotted and broken lines, respectively. 



and three-magnon (l/T^) contributions are also 



APPENDIX A: SPIN- WAVE HAMILTONIAN 



Performing the Fourier and then Bogoliubov transformations as N 1 / 2 J2 n e lfc (™ 1//4 ^ajj = a\.ch8k — k sh8 k and 

Ar_1/2 E„ e " ifc( ™ +1/4)6 ™ = 4 ch(? * - a k shd kl where we abbreviate coshflfc and sinh6>fc as chO k and sh6 kl respectively, 
the Hamiltonian is represented as 



Hi = -2JN[2Vslr - (S + s)A] - [g s S - g s s - (g s ~ g s )A]^HN 



(Al) 



Ho = -2JN\r 2 + A 



\r 2 + a 2 - ( 



S is 



rA 



— ^2 5(fa - k 2 - k 3 + fa) V (1 \k 1 ,k2,k3,k 4 )ai i ai i a k2 a k3 + V {2) (fa,k 2 ,k 3 ,k i )0\ i 0l i k2 k3 



AN 



k1.---.k4 



+V r (3) (^1,^2,^3,^4)0^ 0l 2 a k3 ki - V (4) (fa ,k 2 ,k 3 , fa) (a^ a k2 a ka kl + H.c. ) 
-V Q (5 \fa,k 2 , k 3l k 4 ,)(0l i k2 k3 a ki + H.c.) + V i6) (fa , k 2 , k 3 , fa){0 kl a k2 k3 a kll + H.< 



where 



k 



2N 



W f(fc) = (S + S + 

= w l (k)±(S-s)T 



2J 

(9s + 9s)vbH 



ch2(9 fc - 2VS'scos-sh26' fe ±{S- s) T 



2J 



7l (fc) = 2%/S'scos-ch26> fc 



2,1 

S + s 



(gs - g s )m H 



2,7 



S . /£ 

s v s 



11 ^~ + ^' 



70 (fe) = 
l/ (1) (fc 1 ,A :2 ,fc 3 ,fc4) = ( 

■t 



r-2A 



ch26» fe + 



S ls_ 
s v S 



sh26» fe , 

A-2r 



cos— sh20i, ± 
2 



cos— ch2# fc 



sh26» fr 



fc 4 - fco 



! ^3 - fa 

cos z h cos — - — 



(ch^sh^ch^sh^ + shfl^ch^sh^clu^ 



+ cos- 



fa - k 3 



k 2 - fa 



cos — - — J (ch6>fc 1 ch6>fc 2 sh0 fe3 sh6>fc 4 + sh6 kl sh6 k2 ch6 k3 ch9 ki 
fa 



(A2) 



(A3) 



(A4) 



(A5) 



sh6> fe3 sh6» fc4 (cosych6> fel sh6> fe2 + cosysh0 fcl ch0fc 2 ) + sh6> fcl sh6» fe2 ^cosych0 fe3 sh0fc 4 + cosysh^ch^^ 



6 



chflfe, ch0. 



fe4 ( C0S y 



fc k 
1 sh0 fcl ch0 fc2 + cosy-ch0 fcl sh0 fc . 



) + ch0 fel ch0 fe2 (c 



2 



ch0 fcl ch0 fc2 ( cos— sh0 fc3 ch0 fc4 + cos— ch0fe 3 sh0fc 4 



&4 

2~ 



V (2) (fci,A:2,fc3,fc4) = (cos 



k 4 - k 2 



k 3 - ki 



cos ) (ch0 fcl sh0fe 2 ch0fc3sh0fc 4 + sh0 fcl ch0 fc2 sh0fc 3 ch0fc 4 



/ fc 4 - fc 3 — fei 
+ (^cos h cos 



2 2 

(ch0 fcl ch0 fc2 sh0 fe3 sh0 fe4 + sh^ sh0fe 2 ch0fc 3 ch0fc 4 



ch0 fc3 ch0 fe4 (cosysh0 fcl ch0 fe2 + cosydi0 fcl sh0 fc2 ) + ch0 fel ch0 fc , 
sh0 fc3 sh0 fe4 (cosych0 A . 1 sh0 fc2 + cosysh0 fcl ch0 fc2 ) + sh0 fcl sh0 fc2 



fc 4 - fc 2 
cos h cos 



&3 — 



(ch0 fcl ch0 fc2 ch0 fe3 ch0* 



/ fc 4 - k 3 k 2 -k 1 
+ cos — h COS 



2 2 

(ch0 fcl sh0 fc2 sh0 fe3 ch0 fc4 + sh0 fel ch0 fc2 ch0 fe3 



y*— " sh0 fe3 ch0 fc4 ^cosych0 fcl ch0 fe2 + cosysh0 fel sh0 fe2 ) + sh0 fcl ch0 fc2 
ch0 fe3 sh0 fc4 ^cosysh0 fcl sh0 fc2 + cosych0 fcl ch0 fe2 ) + ch0 fcl sh0 /C2 



V r 



Q 4: \k 1 ,k2,k 3 ,k i ) = ( 

t 



fc 4 - fc 3 
+ cos — 1- COS 



fc 4 - k 2 

ct)S h COS 

k 2 - fci 



k : 
1 

fc 3 - fej 



(ch0 fcl sh0 fc2 ch0 fe3 ch0 fe 



(ch0 fcl ch0 fc2 sh0 fe3 ch0 fc4 + sh0 fcl sh0 fe2 ch0 fc . 



ch0fc 3 sh0 fc4 ^cosysh0 i;i ch0 fe2 + cosy ch0 fcl sh0 fc2 ) + ch0 fcl ch0fc 
sh0 fc3 ch0 fc4 ^cosych0 fcl sh0 fe2 + cosysh0 /ci ch0 /i;2 ^ + sh0 fel sh0 fe 



V^\k u k 2l k 3 ,k 4 ) = ( 

+( 



; 4 — k 2 k 3 — ki 
cos h cos 



fc 4 - k 3 
cos h cos 



2 

fc 2 - fci 



(ch0 fcl sh0 fc2 ch0 fe3 ch0 fe4 



(ch0 fel ch0 fe2 sh0 fc3 ch0 fe4 + sh0 fel sh0 fe2 ch0 fc3 



sh0 fe3 ch0 fc4 ^cosych0 fcl sh0 fe2 + cosysh0 fcl ch0 fc2 ) + sh0 fel sh0 fc2 

ch0fc 3 sh0fe 4 (cosysh0 fcl ch0 fe2 + cosycb^sb^fc.^ +ch0fc 1 ch0fc ! 
ki - k 3 t k 2 - fci 



cos ] (ch0/ ci sh0 fc2 sh0 fc3 ch0 fe 



^0 (6) (fcl^2,fc 3 ,fc4) = 2(cOS- ^ 

ch0 fc3 sh0fe 4 (^cosysh0 fcl sh0 fe2 + cosych0 fcl ch0 fc2 ) +ch0 fcl ch0 fc 
fci . „ , „ k 2 



sh0 fe3 ch0 fc4 (cosych0 fcl ch0 fe2 + cos Ysh0 /£l sh0 fe2 ^ + sh0 fcl ch0 fc2 



cosysh0 A:3 ch0 fc4 + cosyd^sl^ 

fc fc,i 
cosych0 fc3 sh0 fc4 + cosysh0 fc3 ch0 fc4 

- sh0 fel sh0 fc2 sh0 fe3 sh0 fe4 ) 
sh0 fc4 ) 

cosych0 fe3 ch0 fc4 + cosysh0 fc3 sh0 fe4 
cosysh0 fe3 sh0 /C4 + cosych0fc 3 ch0 fc4 

- sh0 fcl ch0 fc2 sh0 fc3 sh0 fe4 ) 
sh0 fc4 ) 

cosysh0fc 3 sh0 fc4 + cosych0 fe3 ch0 fc4 
cosych0 fc3 ch0 fc4 + cosysh0 fe3 sh0 fc4 

- Sh0 fel ch0 fc2 sh0 fe3 sh0 fe4 ) 
sh0 fe4 ) 

cosych0 fc3 ch0 fc4 + cosysh0fc 3 sh0 fel 

COSySh0fe 3 sh0fe 4 + COSyCh0fc 3 ch0fc 4 



sh0 fcl ch0 fe2 ch0 fe3 sh0 fc4 



COSySh0fe 3 sh0 fc4 + cosych0 fe3 ch0 fc 
cosych0 fe3 ch0 fe4 + cosysh0 fc3 sh0fc. 



(A6) 



0fc is determined so as to diagonalize the linear spin-wave Hamiltonian Tii, that is, to satisfy 71 (fc) = 0. at (f3 k ) 
creates spin waves of ferromagnetic (antiferromagnetic) aspecti^S Minimizing the free energy under the effective 
condition of zero staggered magnetization,-^ W e obtain the optimum distribution functions as 



(a\a k ) = n k , ((i\(i k ) =n+; n% 



1 



e [Ju1(k)-nch28 k ]/k B T _ ±> 



where fi is determined through 



^ ^ S + S + (g s -g a )» B H/2J_ % = {g + s)Ne _ Ja - {0)/knT _ 



k a 



(A7) 



(A8) 
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APPENDIX B: MAGNON SERIES EXPANSION OF THE RELAXATION RATE 



Considering the significant difference between the nuclear and electronic energy scales and assuming the Fourier 
components of the coupling constants to have little momentum dependence^ as ^ n & kn A T n = A T k ~ A T and 
J2 n e lkn Bn = Bl ~ B T (t = — , z), we obtain the Raman and three-magnon relaxation rates as 



Kjn — 1^1^ ( rfc 2 » fc i) "fa + !) 

J l fcl CT=±T=± 

1 {gmU^B-) 2 



dwf(fc) 



dfc 



(Bl) 



fc=fcj 



n(3) 



7">V . 

J l fc 1 ,fe 2 <r=±T=± 

f,CT(J<7\2 -<7 — (7 /- 



E E E 2r w (rC,fe,fa) 2 <nI 2 K r .+l) 



l&TiSJN 2 

+W° s %k l M,Tkf°) 2 nl l nl 2 {nl r « + 1) 



dwf(fc) 



dfc 



dwf(fc) 



dfc 



(B2) 



with a = —(j, where fc 2 a and fcj' 7 CT arc determined through — Lo k °o +Tiuj^/ J = and u) ki 4 
0, respectively, and 

A z A z 
W (fci,fc 2 ) = — ch^ fel ch^ 2 -sh0 fel sh0fc 2 , W ++ {k 1 ,k 2 ) = — sh0 fel sh0 fe2 - ch0 fel ch(9 fc2 



(B3) 



A' 



W (fci, fc 2 , fc 3 ) = — ch6 kl ch9 k2 ch9 k3 



TF+++(fci,fc 2 ,fc 3 ) = —sh9 kl sh9 k2 sh6 k3 



■ (^—ch9 k3+k2 ^ kl - -y/^shflfe. 



+k 2 — fci 



^sh^sh^ ~ ^^^a + fc.-fcl) 

s (fc 3 + fc 2 — fci) — Tiwn/J 

^ (4) (fc 1 ,fc 2 ,fc 3 ,fc 3 + fc 2 -fci) ^- ^ 

^{h + k^k^ + n^/j {^ sM ^-^ - V - ■ ' ' . ) ^ 

S,„ , „ ,„ ^ (5) (fci,fc 2 ,fc 3 ,fc 3 + fc 2 -fci) 

(fc 3 + fc 2 — fci ) — Ti^n / J 
2ffF (2) (fci , fc 2 , fc 3 , fc 3 + fc 2 - fci) / A~ ^ 
Wi"(fca + fca — fci) + 7ma>n/ J \B~ 3 2 1 



— ch6» fel ch6>fe 2 ch0 fe3 



ch^fca+fca-fc 



0- 



A~ 



W+-+QCU fc 2 , fc 3 ) = 2(— sh0 fcl ch0 fca sh0 fc3 



: ^-g3 c h^fc 3 +fc 2 -fci - J —sh9 k3+k2 _ k 



ch9 kl sh9 ko ch9 k 



sv o (h + k 2 - fci,fc 3 ,fc 2 , fci) 
^i(k 3 + fc 2 - fci) — Huj-m/J 



2SV (5) (fc 3 ,fc 3 + fc 2 -fci,fci,fc 2 ) ( A- 



-sh9 



fc3+fc2— fcl 



^-+-(fci,fc 2 ,fc 3 ) = 2(^cM fcl sh0 fe2 ch^ 3 - 

/ A~ is~ \ 

x y-^ch9 k3+k2 - kl - J —sh9 k3+k2 - kl J - 



u>i (k 3 + k 2 — kx) + Huj n / J \B 

I^ch^aW J - 25^)(fc 3 ,fc 3 + fc 2 -fci,fci fc 2 ) 
/ (fc 3 + fc 2 - fci) - huj^/ J 



ch6 



fc3 + fc2— fcl 



^T/ (3) (fci,fc 2 ,fc 3 ,fc 3 + fc 2 -fci)^ 



OJi (fc 3 + fc2 — fcl) + hui-^/J 



(^^p S Wfc3 + fc2-fcl - W ^Ch^fcj+fcj-fc^. 

(B4) 
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